Disorder-driven splitting of the conductance peak at the Dirac point in graphene 
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The electronic properties of a bricklayer model, which shares the same topology as the hexagonal 
lattice of graphene, are investigated numerically. We study the influence of random magnetic-field 
disorder in addition to a strong perpendicular magnetic field. We found a disorder-driven splitting 
of the longitudinal conductance peak within the narrow lowest Landau band near the Dirac point. 
The energy splitting follows a relation which is proportional to the square root of the magnetic 
field and linear in the disorder strength. We calculate the scale invariant peaks of the two-terminal 
conductance and obtain the critical exponents as well as the multifractal properties of the chiral 
and quantum Hall states. We found approximate values v ~ 2.5 for the quantum Hall states, but 
1/ — 0.33 ±0.1 for the divergence of the correlation length of the chiral state at -B = in the presence 
of a strong magnetic field. Within the central n — Landau band, the multifractal properties of both 
the chiral and the split quantum Hall states are the same, showing a parabolic /[a(s)] distribution 
with a(0) = 2.27 ± 0.02. In the absence of the constant magnetic field, the chiral critical state is 
determined by q(0) = 2.14 ± 0.02. 

PACS numbers: 73.23.-b, 71.30.-|-h, 73.22.-f 



I. INTRODUCTION 

The nature of the current carrying states near the 
charge-neutral Dirac point in graphene is of exceptional 
interest and substantial importance for the understand- 
ing of the electrical transport properties in strong mag- 
netic fields. The experimental observation of quantum 
Hall plateausi-^i^iii^i^ with a^y = {2N + 1) 2e^/h, N = 
0, 1, . . . strikingly emphasizes the significance of disorder, 
which influences the single sheet of carbon atoms forming 
the hexagonal lattice of graphene. Due to the two val- 
leys appearing in the band structure, each Landau band 
contributes two times 2e^/h to the Hall conductivity cTxy, 
whereas the factor of 2 accounts for the spin degeneracy. 

In the absence of a magnetic field, simple on-site (di- 
agonal) disorder gives rise to strong Anderson localiza- 
tion, which causes the electrical current to vanish at zero 
temperaturesi^ This is, however, in conflict with exper- 
imental observationsii^i^ii^ demonstrating that either a 
different type of disorder is present in real samples or that 
electron-electron interaction renders the one-particle pic- 
ture obsolete. A disorder type being able to account for 
a finite conductance, is the ripple-disorderi2ii£ which is 
believediiii^ to create similar effects as those originat- 
ing from fluctuations of the hopping terms due to elastic 
strains of the intrinsic curvatures of the graphene sheet. 
It is also well known that the inter-valley scattering de- 
pends crucially on the type of disorderi ^^i^^i^^i^^i^^ 

Recently, a splitting of the conductivity maximum 
within the central Landau band at the Dirac point has 
been observed in high mobility graphene samples for very 
strong magnetic-flux densities B > 20 Ti^i^i^ii^ The mea- 
sured energy splitting AE oc \/B has been suggested^ to 
be due to a lifting of the sub-lattice symmetry caused 
by electron-electron interaction. Also, the effect of 
counter-propagating chiral edge states^^, electron-lattice 
effectsr^ as well as valley ferromagnetism22i2ii^i22ii^ 



have been put forward to account for the observed split- 
ting near the Dirac point. Although the proposed ap- 
proaches contain many interesting physics based on inter- 
action effects, we would like to retain the non-interacting 
particle picture in the present work and investigate the 
influence of a random magnetic field (RMF) which causes 
similar effects as ripple-disorder . i ^ i ^ ^ Recently, the influ- 
ence of real random hopping terms was studied and a 
splitting of the extended state in the n = Landau band 
has been reported.— However, the experimentally ob- 
served extremely narrow Landau band^i^i^ii^ at the Dirac 
point suggests the origin of the Landau-level broadening 
to be ripple-disorder or an equivalent disorder that also 
preserves the chiral symmetry.— 

By using of a microscopic bricklayer model, which is 
topologically equivalen t^^i^^ to a hexagonal lattice, and 
assuming a random magnetic-flux disorder by introduc- 
ing complex phases into the hopping terms of a tight- 
binding Hamiltonian, we found a narrow density-of-states 
peak at the Dirac point and a splitting of the central 
conductance peak similar to what has been observed in 
experiments. We calculate the density of states, the two- 
terminal conductance, and the critical eigenstates, from 
which the respective energy and magnetic-fleld depen- 
dence of the energy splitting, the scaling of the con- 
ductance, and the multifractal properties of the critical 
eigenstates are obtained. The splitting AE increases lin- 
early with the strength of the random flux amplitude 
and shows a ^/B dependence as observed in experiments. 
Besides the split quantum Hall conductance peak, we 
found a central chiral state a,t E = 0. The latter exhibits 
a critical exponent v = 0.33 ± 0.1, which determines 
the divergence of the localization length ^/^q = l^-l 
whereas the remaining two split bands of the central Lan- 
dau level belong to the ordinary quantum Hall symmetry 
class with a critical exponent i' « 2.5. The calculated 
multifractal properties turn out to depend on the static 



2 



J) f (I) • (!) • j) 



-©- 



e- 



-e- 



-©- 



! • O 



e- 



-©- 



® • (D • ffi • ® 



-e- 



-®- 



e- 



-e- 



• ® • 



-e- 



• » » » 



e- 



^//2 < (f)x,z/{h/e) < f /2 so that the mean value {(j)xz) 
is zero and the variance is The parameter / mea- 

sures the strength of the disorder with a maximal value 
/=!• 

The corresponding tight-binding Hamiltonian, with 
spinlcss fcrmionic particle creation (c^) and annihilation 
(c) operators, on the bricklayer is 



(2) 



FIG. 1: (Color online) The two-dimensional bricklayer lattice 
which shares the same topology as the hexagonal lattice of 
graphene. The bi-atomic unit cell is indicated by the dashed 
rectangle. For the calculation of the conductance, the sample 
is connected to two semi-infinite leads (thin lines) with square 
lattice topology, and periodic boundary conditions are applied 
in the vertical direction (top vertical lines are connected to 
sites in the lowest horizontal line). 



magnetic field. For S = 0, we found a{Q) = 2.14 for the 
chiral state at the Dirac point. This value changes to the 
usual quantum Hall result a(0) = 2.27 in the presence of 
a spatially constant magnetic field. 



II. BRICKLAYER MODEL AND TRANSFER 
MATRIX METHOD 

Graphene can be represented by a tight-binding Hamil- 
tonian defined on a two-dimensional honeycomb lattice. 
In our numerical calculations, the honeycomb lattice is 
transformed into a bricklayer lattice as shown in Fig. [TJ 
Each site is connected by three bonds with its nearest 
neighbors and has the same topology as the hexagonal 
lattice<21iS For investigations of the spectral properties, 
the differences in bond length and bond angle do not 
matter. In other cases the length scale has to be properly 
adjusted, e.g., for the plaquette size 2a^ (bricklayer) and 
(3V3/2)a^ (hexagonal lattice), where a is the respective 
nearest-neighbor distance. 

A magnetic field perpendicular to the xz plane creates 
a magnetic flux through each individual plaquette. 



P 

^x,z = -h/e + 4>x 



(1) 



where p and q are integers which arc mutual prime, and 
(f>x,z is the random flux component. The latter is assumed 
to be due to ripples and to the buckling of the non-planar 
carbon monolayer which causes the magnetic flux to fluc- 
tuate from plaquette to plaquette. Since the plaquette 
size on the bricklayer system is twice that of the square 
lattice, the magnetic-flux density is i? = p/q x /i/(e2a^). 
We choose (fix.z to be uniformly distributed according to 



The complex phase factors along the vertical bonds of 
a given plaquette are determined by the respective flux, 
&x.z+2a — Sx,z ~ 2tt^x,z- Wc fix the length scale by a = 1 
and set F = 1, which deflncs the energy scale. Please 
note that the first sum in Eq. ^ contains only the 
non-zero vertical hopping terms as shown in Fig. [T] An 
important property of the random flux model is that its 
disorder docs not destroy the chiral symmetry of the sys- 
tem at the band center E = Qi^^^^l^li^l Thus, the Hamil- 
tonian ^ enables us to study both the chiral and the 
quantum Hall critical regimes within the same model. 

For the calculation of the conductance, the samples are 
attached to two semi- infinite leads, which are constructed 
by a two-dimensional square lattice. No magnetic field 
is considered within the two leads. Periodic boundary 
conditions are assumed in the vertical (a;) direction in 
order to eliminate surface effects and edge states in the 
quantum Hall regime. This requires the number of sites 
in the vertical direction to be even, and a multiple of 2q. 

The two-terminal conductance is calculated by the 
well-known transfer-matrix method^^ 



cosh (ei/2) 



(3) 



where A^ch is the number of open channels, t is the trans- 
mission matrix, and the parametrize its eigenvalues. 
Statistical ensembles of typically iVgtat = 10^ samples 
were collected and the mean value of the conductance was 
calculated for each set of parameters {E, L, f,p/q). The 
transformation of the honeycomb lattice into the brick- 
layer changes the length scale along the propagation di- 
rection, which is due to the rectangular shape of the unit 
cell. For computational efficiency, we still consider Lx L 
bricklayer lattices. We verifled that for different shapes 
the conductance scales like gc ^ Lx/L, at all critical 
points. Hence the conductance results we show are al- 
ways by a factor of two larger. Therefore, it is easy to 
correct for our special aspect ratio when comparing with 
conductance results obtained by other researchers. 
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FIG. 2: (Color online) The density of states for L x L brick- 
layer samples with size L — 128, averaged over 500 disor- 
der realizations, showing three Landau bands closest to the 
Dirac point for energy E > 0. The random flux disorder 
strengths are / = 0.001,0.002,0.005,0.007, and 0.015, and 
p/q = 1/32. The inset shows the broadening of the narrow 
central Landau band together with the chiral state at i? = 
for / = 0.001, 0.002, 0.003, and 0.005. 
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FIG. 3: (Color online) The energy dependence of the disorder 
averaged two-terminal conductance {g{E)) for -E > 0, p/q = 
1/32, and various random flux disorder strengths /. The size 
of the square samples is L = 128. Besides the conductance 
peak of the chiral state, which remains at _E = 0, the shift of 
the split quantum Hall state with disorder is depicted. 



Since the number of eigenvalues is the same for each 
principle band shown in Fig. [2l the two split n = sub- 
bands must contain less eigenvalues than the higher Lan- 
dau bands due to the additional chiral states close to the 
Dirac point E = 0. 



III. RESULTS AND DISCUSSION 

A. Density of states 

The density of states p{E) was calculated by counting 
the eigenvalues obtained from diagonalization of = 
128 X 128 samples with periodic boundary conditions in 
both directions, averaged over 500 disorder realizations. 
Without disorder, p{E) was checked to be identical with 
the result of a true honeycomb lattice. Figure [2] shows 
p{E) for p/q = 1/32 flux quantum and several disorder 
strengths / = 0.001, 0.002, 0.003, and 0.005. Due to the 
symmetry around E ^ 0, the density of states is shown 
only for energies E > 0. The Landau spectrum and the 
broadening of the bands are visible. However, compared 
with the n — 1 and n — 2 Landau bands, the central 
n = one remains very narrow and is hardly discernible. 
Therefore, in the inset of Fig. [2l the broadening of the 
narrow central band is shown in more detail where also 
the chiral peak at £^ = is distinguishable in the mid- 
dle of the broadened Landau band. A real splitting of 
the n = Landau band becomes apparent only for very 
small disorder strength /. The total width at half height 
of the central band (neglecting the narrow chiral peak) 
is proportional to \/p/qf- It should be noted that the 
degeneracy of the n = Landau level found in the Dirac 
model, even in the presence of a random magnetic field^ 
is already lifted in our clean lattice model due to Harper's 
broadeningi^^ In the case B = 0, this lattice effect also 
causes deviations from the linear energy dispersion away 
from the Dirac point. 



B. Conductance peak splitting 

The energy dependence of the disorder averaged two- 
terminal conductance is shown in Fig. [3] for different ran- 
dom magnetic-field disorder strengths /. Again, only the 
positive energies are shown, because the conductance is 
an even function of the energy E. With increasing /, the 
conductance peaks corresponding to the n = Landau 
level, with peak value of about 1.12 e^//i, move away from 
E ^ 0, where the conductance peak of the scale indepen- 
dent chiral state remains fixed at g{E = 0) = 1.27e^//i. 
The disorder induced splitting AE of the n = con- 
ductance peak is plotted in Fig. 0] versus / for various 
magnetic fields with p/q = 1/24, 1/32, 1/64, and 1/128. 
For not too large /, the data points follow the straight 
solid lines which are given by the relation 



AE = 




Therefore, the conductance peak splitting shows a ^/B 
behavior as observed in experiments and increases lin- 
early with the strength / of the random magnetic-flux 
disorder. Neither a splitting nor any shift of the conduc- 
tance peak due to the disorder / investigated could be 
observed for the n = I Landau band and L < 256. 

As known from experiment^ the minimum conduc- 
tivity at the Dirac point depends verifiably on disor- 
der. Therefore, in Fig. [5] we show the disorder depen- 
dence of the critical chiral conductance {g{E = 0, /)) 
for different magnetic fields. The two values 4/7r(e^//i) 
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FIG. 4: (Color online) The splitting of the conducting states 
from the n = Landau band obtained from the two-terminal 
conductance of L = 256 square samples versus random mag- 
netic field disorder strengths /. The magnetic fields are 
p/q = 1/24 (T), p/q = 1/32 (A), p/q = 1/64 (■), and 
p/q — 1/128 (♦), respectively. The straight lines follow 
AE=^/^f. 



FIG. 5: (Color online) The averaged two-terminal conduc- 
tance {g{E = 0, /)) for several magnetic fields (given by the 
ratio p/q) versus random magnetic-field disorder strengths /. 
p/q = 1/64: L = 256 {+), L = 512 (O), L = 896 (■). 
p/q = 1/128: L = 256 (x). B = 0: L ^ 256 (□), L = 512 
(O), and L — 896 (A). The horizontal lines mark the con- 
ductance values 4/n {e^ /K) and 8/7r (e^//i), respectively. 



and %/n{e'^/h) are indicated by solid horizontal lines. 
The conductance turns out to be scale invariant for suf- 
ficiently strong random flux disorder. In the limit of 
/ 1.0, {g{E = 0)) converges to a common value that 
is a little larger than S/tt (e'^/h), which is independent of 
magnetic field. In contrast, a magnetic field dependence 
is seen for small /. In the range 0.001 < / < 0.4, the 
conductance seems to approach a value somewhat smaller 
than 4/7r(e^//i) for all magnetic fields p/q < 1/32 stud- 
ied. For small disorder, e.g., / = 0.1 and i? = 0, the 
conductance scales with the sample size. A convergence 
to a value sa 4/7r (e^//i) as in the finite B case is compat- 
ible with our data. 

In a finite clean system, {g{E = 0, / = 0)) is zero for 
B = but increases with the system size for finite B. 
The former behavior is due to the vanishing density of 
states at E ~ in the absence of a magnetic field. In 
the thermodynamic limit, the result for the conductance 
seems to depend on the order of limits L — > oo and f 0. 
A special case appears in the clean limit if is a multiple 
of three. {g{E = 0, / = 0)) turns out to be 2e^/h due 
to the four eigenstates appearing at = for /c^ = tt/S, 
where = Ntt/Lz and N = 1,2,... . This agrees well 
with the experimentally observed minimal conductivity 
fmin = 4e^//i if the electron spin is taken into account. 



C. Critical quantum Hall regime 

The energetic position of the split conductance peak 
within the n = Landau band as shown in Fig. [3] was 
estimated from the energy and size dependences of the 



conductance peak. For all values of disorder strength /, 
the conductance g{E) exhibits a maximum around the 
critical energy Ec{f). The maximal value of the con- 
ductance is gc ~ 1.12e^//i. If we take into account the 
factor of two due to the special aspect ratio, this value is 
in agreement with the critical conductance gc = 0.6e^/h 
found on the square lattice with spatial correlated diago- 
nal disorder!^ Likewise, it compares well with the results 
gc = 2/Tre'^/h obtained within a self-consistent Born ap- 
proximation by Shon and Ando^ for a graphene sheet 
in the presence of short- or long-range scattering poten- 
tials. Finally, our result for the split conductance peak 
agrees also with the analytical calculations of Ostrovsky 
and co-workers ji^ii^ii^ 

A more detailed analysis of the size and energy de- 
pendences of the conductance shows that g{E) exhibits 
a pronounced maximum around the disorder dependent 
critical energy Edf). This renders the estimation of 
the corresponding critical exponent possible. The peak 
g{E^ /) becomes narrower when the size L of the sys- 
tem increases. Following the scaling theory of localiza- 
tion, the conductance is assumed to be a function of 
only one parameter in the vicinity of the critical point, 
{g) = F[L/S^{E)\, with a diverging correlation length 
S,{E) (X \E — Ec\~'^ ■ Then, we can extract the critical 
value of the conductance and the critical exponent as- 
suming that 

(.g) ^gc- A{L) [E - Ecf . (5) 

and 

A{L) oc L^/''. (6) 
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FIG. 6: (Color online) The energy dependence of the averaged 
conductance {g) for p/q = 1/32 and / = 0.10 calculated for 
the bricklayer model of size L x L. Solid lines are quadratic 
fits, given by Eq. ([5]). A scaling analysis provides us with a 
critical exponent f ~ 2.5 if only the data for L > 384 were 
used. The inset shows a double-logarithmic plot of the size 
dependence of A(L) defined in Eq. ([5]). The solid line is a 
power-law fit of Eq. Q with ~ 2.5. 



Figure [S] shows the size and energy dependence of the 
conductance peak for disorder / = 0.1. We obtained a 
critical exponent w 2.5 when only the data for L > 384 
were used in the scaling analysis. We observed a similar 
behavior for smaller disorder strength / = 0.01. In view 
of the uncertainties, we consider this to be consistent with 
the well-known results v = 2.3 for ordinary quantum Hall 
systems supporting the view of the universality of the 
continuous quantum Hall phase transition. 

The two-terminal conductance for the n = 1 Landau 
band is shown in Fig. [7] versus energy. The peak value 
is about 1.87e^//i. The maximal available size of the 
system (L = 512) is still insufficient for a more accurate 
estimation of the critical parameters. Therefore, we can 
neither give any limiting value nor address the Landau- 
level dependence of the critical conductance. The very 
slow L dependence of the conductance in the critical re- 
gion, indicating a larger localization length than for the 
n = case, makes the scaling analysis and estimation of 
critical exponent impossible. 



D. Critical chiral state 

From the energy and size dependences of the mean 
conductance in the vicinity of the chiral critical point we 
extract the critical conductance, gc ~ limL_,oo(.9(-E' = 
0,L)). Our data for the mean conductance (g) at the 
band center E ^ exhibit a convergence to the size in- 
dependent limit gc, 
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FIG. 7: (Color online) The disorder averaged two-terminal 
conductance versus energy in the n — 1 Landau band for var- 
ious system sizes L = 128, 256, 384, and 512. The magnetic 
field is p/q = 1/32 and the disorder strength / = 0.01. 



As shown in the inset of Fig. El (?c — 1.267 e^//i. 
This perfectly agrees with the theoretically predicted 
value A/ire'^/h = 1.273 e^//i for the spin-degenerate 
situation.— ii^ii^ We mention again that in our case the 
factor of two originates not from the spin but from the 
special aspect ratio used. 

To estimate the conductance fluctuations, we calcu- 
lated the probability distribution p{g) which proved to 
be Gaussian, confirming the presence of mesoscopic fluc- 
tuations also at the chiral critical point. The variance 
var g = {g"^) — {g)"^ varies slightly as a function of B and 
/. For instance, we found var g w 0.23 (e^//i)^ for the 
data shown in Fig.[8]but 0.21 [e^/hf for the 5 = case 
(Fig.©. 

To analyze the critical behavior of the conductance 
in the vicinity of the band center, we calculated {g) for 
square samples of size L = 64, 128, 192, 256, and 384, 
within the narrow energy interval 10~^ ^ E < 10~^. 
Our data confirm that, after finite-size correction ([7]), 
the mean conductance is a function of one parameter 
only, {g{E,L)) = G{EL^/^). Following scahng theory, 
v is the critical exponent that governs the divergence of 
the correlation length, £,{E) oc A more detailed 

analysis shows that 



G{x) ^ gc^ cq ln(l + cix + C2x'^), 



(8) 



{g{E = Q,L))^gc + b/L. 



(7) 



with CO = 0.05, ci = -0.02944, and C2 = 0.02134. Scaling 
analysis provides us with the critical exponent v = 0.33± 
0.1. Although the best scaling is obtained for v = 0.33 
(as shown in Fig. [8]), a reasonable scaling is also possible 
for values of the critical exponent in the range 0.25 < 
V < 0.45. Again, only numerical data for much larger 
systems would reduce this inaccuracy. 
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FIG. 8: (Color online) Scaling of the averaged chiral con- 
ductance in the vicinity of the critical Dirac point E — 0. 
The magnetic field is p/q = 1/32 and the disorder strength 
/ = 0.01. The data show scaling with a critical exponent 
u = 0.33 ± 0.1. The solid line is the fitted function p. The 
inset shows the finite-size correction g{E — 0, L)/{e^ /h) — 
1.267 + 3.178/L, and var g ^ 0.23 (e^/h)^. 
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FIG. 9: (Color online) Attempt to scale the averaged conduc- 
tance data close to the chiral critical point for a bricklayer 
system with zero homogeneous magnetic field and random- 
ness / — 1.0. Bricklayer square samples of size L — 128, 
192, 256, 384, and 512 were considered. We obtain a crit- 
ical exponent u = 0.42. The inset shows the size depen- 
dence of g{E = 0) for 128 < L < 960. The solid line is 
g{L)/{e^/h) = 2.62 - 11.48/L, and var g w 0.21 {e^/hf. 



We have carefully checked that our estimate of v is 
insensitive to the range of the energy interval selected. 
Also, we tried to compensate for the observed divergence 
of the critical density of states near E = 0, which, how- 
ever, changed our estimates only marginally. We mention 
that our result does not satisfy the Harris criterion^ 
which states that dz^ — 2 > 0. There are similar results 
also for other model o^°i'^"'^'^^i^^''^'^i'^^ with chiral critical ex- 
ponents 1/ < 1. 

In the absence of a constant magnetic field, we ob- 
served that the value of the conductance at = de- 
pends considerably on the size of the system, g{L) = 
2.62 — 11.48/L. Even worse, the scaling behavior is ful- 
filled only in a very narrow interval of conductance val- 
ues. From Fig. [H] we see that the width of the scaling 
interval is comparable with the uncertainty in the con- 
ductance due to the finite-size effects. It is therefore no 
surprise that the estimation of the critical exponent is 
much more difficult. Taking into account these restric- 
tions, we conclude that the estimated value of the critical 
exponent i/ ~ 0.42 is still in reasonable agreement with 
the previously obtained values v = 0.35 and ly = 0.42 for 
a normal square lattice^ii^ 



E. Multifractal eigenstates 

The critical eigenfunctions ipEif) of bricklayer samples 
with sizes up to = 384 x 384 were obtained numerically 
using a Lanczos algorithm. For eigenvalues close to E' = 
0, the chiral eigenstates show the expected sub-lattice 



polarization, since for a clean system the wave functions 
are non-zero only on one of the sub-lattices,^>i^ This 
sublattice polarization is still observable in the presence 
of both the fluctuating random flux and the constant 
perpendicular magnetic field. We also notice some weak 
quasi-one-dimcnsional structures making the appearance 
of |'(/'E=o(7')P look slightly anisotropic. 

A similar behavior is found for the n ~ critical quan- 
tum Hall states, as well. This is shown in Fig. [TOl where 
the squared amplitude of a characteristic eigenstate with 
energy E = 0.016815 is displayed for a L = 320 sample 
with p/q = 1/32 and RMF strength / = 0.1. A closer in- 
spection reveals an approximate sub-lattice polarization, 
also in this case. This effect appears to be stronger than 
in the diagonal disordered situation4^ 

The multifractal analysis was carried out as usual by 
utilizing the well-known box counting method where the 
scaling of a "box-probability" is calculated, P{s, A) = 

Ef ^"(E.60.(/) \^E{r)n - A-(^), from which the gen- 
eralized fractal dimensions D(s) ~ r(s)/(s — 1) or, by 
a Legendre transform, the so-called /[a(s)] distribution 
can be derived42ii^ Here, fti{l) is the ith box of size 
I ~ XL from which the sth moment of the modulus of 
the normalized eigenstate t/jEir) is taken. 

In Fig. [11] we show the /[a(s)] distribution of two 
d ~ 2 chiral eigenstates at i? « in comparison with 
the parabolic approximation, /[q;(s)] = d — [oi{s) — 
Q!(0)]^/{4[q;(0) — d]}, which for finite systems is nor- 
mally valid only for small |s|. We found a value aifi) — 
2.14 ± 0.02 when averaged over several critical chiral 
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FIG. 10: (Color online) The probability density of a charac- 
teristic critical quantum Hall eigenstate of the n — Landau 
band at E = 0.016815. The system size is L = 320, con- 
stant magnetic field is p/q = 1/32, and random magnetic 
field strength is / = 0.1. 

states in the case B = 0. This value turns out to be 
close to the one obtained previously for a square lat- 
tice with correlated random magnetic field4^ For finite 
B, we found an averaged a{Q) = 2.27 ± 0.02 which is 
the same as the a{0) value of the n ~ quantum Hall 
states. This result also agrees with a(0) values pub- 
lished for various quantum Hall models in the range 
2.26 < a(0) < 2.29..^'^'^'^ For the n = 1 quantum 
Hall state the achievable system sizes were not sufficient 
for the calculation of a conclusive multifractal spectrum. 

We finally mention that we also checked the E ~ 
2.9725 critical cigcnstates within the Landau band closest 
to the tight-binding band edge in the presence of diag- 
onal disorder. These outer states do not belong to the 
Dirac fermion sequence that can be found only in the 
energy range —1.0 < E < 1.0. Again, we obtained the 
ordinary quantum Hall multifractality with a(0) ~ 2.29 
and a critical conductance gc — l.Oe^/h (including the 
aspect ratio factor of 2). 



IV. CONCLUSIONS 

We investigated numerically a two-dimensional brick- 
layer lattice model, which shares the same topology as 
graphenc's hexagonal lattice, in the presence of both a 
homogeneous and a spatially fluctuating random mag- 
netic field. We calculated the one-particle density of 
states, the two-terminal conductance, and the multifrac- 
tal properties of critical eigenstates. Within the very 
narrow n = Landau band, we found a splitting of 
the two-terminal conductance peak into three subpcaks. 
A central chiral conductance peak located at the Dirac 




FIG. 11: (Color online) The f[a{s)] distributions for two chi- 
ral eigenstates with energy close to i? = 0. The parabolas are 
determined by the values a(0) — 2.14 in the case of B = 0, 
/ = 1.0 (x), and a(0) = 2.265 for p/q = 1/32, / = 0.5 (+), 
respectively. In both cases, the system size is L = 320. The 
data points correspond to s = ±4.0, ±3.0, ±2.5, ±2.0, ±1.5, 
±1.0, ±0.5, and 0.0. 

point E ^ Q with a scale independent value A/ire'^/h < 
gc{E = 0,/) < S/ne'^/h (including the aspect ratio fac- 
tor of 2), depending on the strength of the random mag- 
netic field /. The symmetric splitting of the two other 
peaks AE{p/q, f) increases with the square root of the 
applied perpendicular magnetic field expressed by the ra- 
tional number p/q times a flux quantum h/e per pla- 
quette, and linearly with the amplitude of the random 
magnetic- flux disorder /. The scale independent conduc- 
tance peak value of these critical quantum Hall states is 
gc ~ 1.12 e^/h (including the aspect ratio factor of 2), 
independent of disorder strength and applied magnetic 
field. The splitting of the n = Q conductance peak al- 
lows for a Hall plateau with both a^y = and axx = 0, 
except at E = 0, where the longitudinal conductance is 
finite due to the chiral critical state. 

A similar scenario, a splitting of the Landau band and 
an additional central critical state with a IndiJp singu- 
larity in the density of state was previously reported for 
a two state Landau model by Minakuchi and Hikami^. 
They found divergences of the localization length with an 
exponent v = 0.26 for the central state, which is not far 
away from our result for the chiral critical state, and for 
the split Landau band ly ~ 3.1, which is somewhat larger 
but still believed by these authors to be compatible with 
the conventional quantum Hall universality class. Pre- 
sumably, the large value has to be attributed to finite-size 
effects. 

Our analysis of the critical eigenstates revealed the 
sublattice polarization, known from the clean system at 
E = 0, to exist at least approximately also in the random 
magnetic-flux disordered system for energies in the vicin- 
ity of the Dirac point. Within the achieved uncertainty, 
the multifractal properties of the critical eigenstates in 
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the n = Landau band appeared to be the same. We 
found the quantum Hall critical states parabolic /[a(s)] 
distribution, which is determined by an a(0) = 2.27±0.02 
value, to be in accordance with the multifractal proper- 
ties of the chiral critical state close to the Dirac point. 
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